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ABSTRACT 

Very massive high redshift clusters can be used to constrain and test the ACDM 
model. Taking into account the observational constraints of Jee et al. (2009) we have 
calculated the probability for the most massive cluster to be found in the range (5.2 — 
7.6) x 10 14 Mq, between redshifts 1.4 ^ z ^ 2.2, with a sky area of 11 deg 2 and under 
non-Gaussian initial conditions. Clusters constrain the non-Gaussianity on smaller 
scales than current cosmic microwave background or halo bias data and so can be 
used to test for running of the non-Gaussianity parameter /nl- 
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1 INTRODUCTION 

In this article we use galaxy cluster observations to 
place constraints on primordial non-Gaussinaity. The weak- 
lensing analysis of the galaxy cluster XMMU J2235.3 — 2557 
(hereafter XMM2235) presented by Jee et al. (2009) indi- 
cates the existence of a massive cluster which has a mass 
of M 324 = (6.4 ± 1.2) x 10 14 A/ Q at a redshift z » 1.4 * 
The probability of observing a cluster with mass equal to or 
above the observed one is of the order of 5 x 10~ 3 for the 
surveyed area (11 deg 2 ) and redshift range 1.4 ^ z ^ 2.2. 

The expected number of high redshift clusters is derived 
from the distribution of the primordial matter density fluc- 
tuations. For canonical single field inflation models the pri- 
mordial matter density fluctuations are effectively Gaussian 
distributed. However, various alternatives have been pro- 
posed in which different amounts of non-Gaussianity are al- 
lowed (see for example Komatsu et al. (2009) for a review). 
Deviations from Gaussianity are often parametrized by the 
dimensionless non-linear coupling parameter /nl (Komatsu 
& Spergel 2001) through the 'local' relationship 

*(x) = 0(x) + f NL (0(x) 2 - (0(x) 2 )), (1) 

where 4> denotes the Bardeen gravitational potential and <f> 
represents an isotropic Gaussian random field. Several works 
have calculated the effect of /nl 7^ on the predicted cluster 

* The subscript indicates that the halo is denned as a sphere 
with a density which is 324 times the mean matter density of 
the Universe. Also, all errors will be given at the one standard 
deviation confidence intervals unless stated otherwise. 



number counts (Matarrese et al. 2000, Lo Verde et al. 2008, 
Grossi et al. 2009, Desjacques, Seljak & Iliev 2009, Des- 
jacques & Seljak 2010). It has been pointed out by Jimenez 
& Verde (2009) and Sartoris et al. (2010) that for cluster 
XMM2235 not to be a rare event, values of /nl above the 
current limits set by CMB observations might be required. 
Sartoris et al. (2010) included a discussion on the degener- 
acy between /nl and as- Conclusions from both works are 
based on the expected number of clusters with masses equal 
to and above 5 x 10 14 Mq at the observed redshift. 

Holz & Perlmutter (2010) proposed evaluating the prob- 
ability of the most massive observed cluster in a survey. As 
the most massive cluster will be one of the easiest to find 
in the survey, it will be less sensitive to selection effects. 
They find that XMM2235 is near the 3 sigma error con- 
tours in maximum mass and redshift for a 11 deg 2 survey. 
Note that assuming XMM2235 is the most massive cluster 
in the survey is conservative in that if it turns out there 
is a more massive one in that survey the probability would 
be even lower. In this article, we extend Holz & Perlmut- 
ter's approach by evaluating how the probability of the most 
massive cluster is affected by non-Gaussianity. 

The rest of this article is organized as follows. The the- 
oretical background for calculating the expected number of 
counts in the presence of non-Gaussianity is presented in 
Section 1. A discussion on how the probability distributions 
are built and what are the results obtained for the observed 
range of masses and redshifts is given in Section 2. The con- 
clusions are given in Section 3. 
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2 NUMBER COUNTS IN THE PRESENCE OF 
NON-GAUSSIANITY 

Based on the Press- Schechter formalism, Matarrese et al. 
(2000) and LoVerde et al. (2008) computed analytical ex- 
pressions of the mass function dn^o/dM in the presence 
of non-Gaussian initial conditions. These expressions have 
been tested against simulations (Grossi et al. 2009) for 
masses up to ~ 10 15 Mq and /nl ~ 100. They indicate 
that the Matarrese et al. (2000) approach produces a better 
fit at high masses/larger /nl (also see G. D'Amico et al. 
2010, Maggiore & Riotto 2009). For this reason we use the 
Matarrese et al. (2000) approach for modeling the effects of 
non-Gaussianity, although we use it somewhat beyond the 
ranges it has been tested for by N-body simulations. 

One can write the non-Gaussian mass function in terms 
of the Gaussian one dn/dM multiplied by a correction factor 



dn^ G (M, z, /nl) _ dn(M,z) 



dM 



p m dlna(M,z) 



dM 



Rng{M, z, /nl)) 



M 



dM 



-f(a(M,z))R NG (M,z,f NL ), 



where p m = fl m p c , a 2 {M,z) is the variance of the 
smoothed density fluctuations and f(a(M,z)) is the mass 
function. We assume non-Gaussianity deviations of the 'lo- 
cal' form and adopt the CMB convention in the definition 
of /nl- The suggested barrier factor of 0.75 is introduced in 
the computation of the critical collapse density to account 
for deviations from spherical collapse (Grossi et al. 2009). 
Rng increases with the skewness of the density distribution 
which increase for larger values of /nl ■ The above expression 
depends as well on the transfer function for which we use the 
fit by Bardeen et al. (1995) with the correction introduced 
by Sugiyama (1995). We checked that the predicted number 
of clusters is not sufficiently different when the more accu- 
rate Eisenstein & Hu (1998) formula is used for our results 
to be significantly affected. 

We are ultimately interested in calculating the prob- 
ability density function of the most massive cluster within 
certain mass and redshift ranges. Such a calculation is based 
on the predicted average number of clusters. The average 
number of clusters per redshift interval dN/dz in the mass 
range [M min ,M max ] is given by 

"M max 



dN(z, / N 
dz 



f. 



sky " 



dV(z) 



dz 



dM 



dn NG (M, z, /nl) 
dM 



(2) 

where / s k y represents the fraction of sky observed, dV(z)/dz 
is the volume element given by 



ay_ 

dz 



47T 

Wz) 



dz' 



H(z 



and H(z) is the Hubble parameter. 



3 PROBABILITY OF THE MOST MASSIVE 
CLUSTER UNDER NON-GAUSSIAN 
INITIAL CONDITIONS 

We were interested in answering the following question: 
What is the probability of galaxy cluster XMM2235 be- 
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Figure 1. Variation of the probability of the most massive clus- 
ter, for different mass bins, as a function of non-Gaussianity (/nl) 
and erg . The blue (dashed), green (solid) and red (dash-dot) curves 
are different mass bins and the thick and thin curves are different 
a-g values. 



ing the most massive cluster within the surveyed volume? 
Jee et al. (2009) estimated the mass of this cluster to be 
M 324 = (6.4 ± 1.2) x 1O 14 M . We use the ACDM Jenk- 
ins et al. (2001) mass function fitted for the 324 overden- 
sity halo definition. We discuss the effect of mass function 
choice in Section 4. The redshift interval of the survey was 
1.4 z ^ 2.2 and the area 11 deg 2 (f a k y )- We took as 
our fiducial cosmology, the best fit WMAP7 cosmological 
parameters (Larson et al. 2010): h — 0.710; Q m = 0.266; 
Q. b h 2 = 0.02258; n s = 0.963; erg = 0.801. We had three mass 
bins: M324 < 5.2 x 10 14 , 5.2 x 10 14 < M 32 4 «S 7.6 x 10 14 , 
and M324 > 7.6 x 10 14 . By making our bin size of order or 
larger than the error bars of the observed mass estimate, 
we approximately account for the effects of Eddington bias 
(Mortonson et al. 2010). We assumed the number of clusters 
in each mass bin was Poisson distributed with the expected 
number of clusters given by Eq. 2 integrated over the red- 
shift range of the survey. We then simulated a large number 
of surveys and evaluated how many times each mass bin 
contained the most massive cluster of the survey. The prob- 
ability of the most massive cluster being in each mass bin 
is then worked out by dividing by the total number of sur- 
veys simulated. In Fig. 1 we show how the probability of 
the largest cluster varies with non-Gausianity and erg. As 
would be expected, large non-Gaussianity and large erg lead 
to larger most massive clusters being more probable (Des- 
jacques et al. 2009, Desjacques & Seljak 2010, Sartoris et 
al. 2010). This means that too much non-Gaussinity or too 
high a value of erg is also disfavored as then the most mas- 
sive cluster is more likely to be even larger than XMM2235. 
This is another aspect that distinguishes our approach from 
that adopted by Jee et al. (2009), Jimenez and Verde (2009) 
and Sartoris et al. (2010). They looked at the probability 
of getting a more massive cluster than XMM2235. In their 
approach, increasing the non-Gaussianity or erg will always 
lead to an improvement in the probability. 

We found that our results were very insensitive to the 
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other cosmological parameters in the range of variation al- 
lowed by WMAP7 and so we kept them fixed at their 
WMAP7 maximum likelihood values, except for as. How- 
ever, it should be noted that the mass estimate of (Jee et al. 
2009) depended on the assumed concentration of the cluster. 
They used a concentration parameter (estimated from the 
N-body simulations of Gao et al. (2008)) of c = 3.20 ± 0.75. 
The concentration parameter can be related to the halo for- 
mation redshift (zf) which for a halo of mass (M), observed 
at redshift (z ), can be defined as the redshift that a fraction 
F of the halo was in halos of mass greater than {fM) where 
/ and F are positive constants less than one (Lacey and 
Cole 1993). Using the extended Press-Schechter formalism 
(Bond et al. 1991), this condition can be expressed as 



erfc 



F 



(3) 



where u)(z) is the growth function modified threshold over- 
density of spherical collapse at redshift z and can be approx- 
imated by (Lacey & Cole 1993; Nakamura & Suto 1997) 

(12 x vr) 2 / 3 

ui(z)D(z) = 3 x i ^ — x (1 - 0.0123 x log 10 [l + x 3 ] ) 



where x - 



20 

(ii^-'-i) 1 / 3 
1+3 



(4) 

and D(z) is the linear theory growth 



function normalized to equal one at redshift zero. The char- 
acteristic halo density is related to the concentration param- 
eter by 

„3 



200 



log[l + c] - c/(l + c) 



(5) 



Navarro et al. (1996) proposed that 5 C was proportional to 
the mean density of the Universe at the time of the halo 
formation. This model was re-calibrated by Gao et al. (2008) 
who found a good fit (for flat cosmologies) was given by 
setting / = 0.01 and F = 0.1 in Eq. 3 and setting 



5 C = 600 



n m (i + z f ) 3 



(i - n m ) + n m (i + zof 



(6) 



The primary dependence of c on the cosmological param- 
eters is due to the dependence of a in Eq. 3 on as- For 
our parameter choices we find that « 2 and so the 
WMAP7+BAO+H0 prior we use for cr 8 (standard devia- 
tion of 0.024) leads to an uncertainty of in c of about 0.048. 
Additionally, the Gao et al. (2008) simulations are based on 
a og = 0.9 cosmology, while the WMAP7+BAO+H0 prior 
we use has a mean of 0.809. This leads to an error of 0.2 
in c. However, this is negligible compared to the N-body 
fitting error for c of 0.75 used by Jee et al. (2009) to es- 
timate the mass error of XMM2235. So, within the tight 
constraints of WMAP7+BAO+H0 we can neglect the de- 
pendence of c on the cosmological parameters. Hence our 
likelihood for XMM2235 alone will not be strictly correct 
outside the ranges of the WMAP7+BAO+H0 prior that we 
use. However, this will not affect our final constraints as they 
will be based on the posterior probability which includes the 
WMAP7+BAO+H0 prior. 

Another systematic that should be considered is the 
effect of contributions from the triaxial shapes of cluster- 
sized halos and uncorrelated large-scale matter projections 
along the line-of-sight. This generally leads to an additional 
mass error of about 20% (Becker and Kravtsov (2010)). To 
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Figure 2. Likelihood (P) that XMM2235 is the most massive 
cluster as a function of non-Gaussianity (/nl) an( i CT 8- 



account for this we added this error in quadrature to the 
error quoted by Jee et al. (2009). This widened our central 
bin width from 2.4 to 3.5 x 10 14 Mq. We found this increased 
bin width had a negligible impact on our results. 

We generated a grid of probabilities in ag and /nl 
for the XMM2235 being the maximum mass cluster. We 
then interpolated this grid to get our likelihood func- 
tion, p(XMM2235 most massive | /nl, erg), which is plotted in 
Fig. 2. As can be seen there is a strong degeneracy between 
/nl and erg as expected from Fig. 1 and the above discussion. 

The WMAP7 data constrains the local non-Gaussianity 
parameter to —10 < /nl < 74 at the 95% CL (Ko- 
matsu et al. 2010). But this non-Gaussianity constraint is 
at wave-numbers of about 0.04 Mpc -1 (Komatsu & Spergel 
2010, Lo Verde et al. 2008, Sefusatti et al. 2009). Non- 
Gaussianity has also been constrained to similar levels and 
scales as WMAP7 by halo bias (Slosar et al. 2008). The 
non-Gaussianity enhancement, Rng depends on /nl only 
through its effect on the skewness of the smoothed density 
field (Matarrese et al. 2000). The skewness depends on an 
integral over two wave numbers, k and k' , and the angle 
between them. We integrate over the angle and plot the 
kernel in Fig. 3 for the density field smoothed on scales cor- 
responding to the mass of XMM2235. As can be seen the 
mean of the kernel is at around k w 0.1 Mpc -1 . Also the 
amplitude of the kernel remain non-negligable up to about 
k « 0.2 Mpc -1 . Therefore, the XMM2235 measurement 
probes smaller length scales than the WMAP7 measurement 
of /nl- There are multi- field models which produce a scale 
dependent /nl with the spectral index generally being first 
order in the slow roll parameters, but larger values are pos- 
sible (see for example, Byrnes et al. 2009). Also, there are 
models which produce a step-function like change in /nl 
(Riotto and Sloth, 2010), so /nl can be small on WMAP7 
scales but large on cluster scales. 

The WMAP7 data (Larson et al 2010) constrains as 
and the other 5 cosmological parameters (Q m , fit., h, n, r) of 
the fiat ACDM model to very high precision. Even tighter 
constraints are obtained when WMAP7 is combined with 
baryon oscillation (BAO) and Hubble expansion data (HO) 
(Komatsu et al. 2010). The marginalized distributions of the 
6 parameters are very close to Gaussian with small degen- 
eracies between the parameters. Also, as the constraint on 
as is much tighter than the one obtainable from XMM2235, 
we can approximate the WMAP7+BAO+H0 as a Gaussian 
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Figure 3. The kernel of the integral for the skcwness of the linear 
density field smoothed on scales of 16 Mpc -1 . The angle between 
the wave-numbers k and k' has been integrated over. 
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Figure 4. The 95% posterior probability function contours (solid 
blue) for WMAP7+BAO+H0 and XMM2235. The 95% prior 
probability function contours from WMAP7+BAO+H0 are plot- 
ted as red dashed lines. The /nl here is on scales much smaller 
than those constrained by WMAP7 or halo bias measurements. 



prior on as with mean of 0.809 and a standard deviation of 
0.024. The posterior for /nl and as is obtained by multi- 
plying the likelihood and prior. The 95% confidence interval 
contours are plotted in Fig. 4. As can be seen by comparing 
Fig. 2 with Fig. 4, there is some tension between the higher 
values of as preferred by XMM2235 and those favoured by 
WMAP7 even when a non-zero /nl is allowed. It can also 
be seen that WMAP7 completely breaks the degeneracy be- 
tween as and /nl- The marginalized posterior probability 
density function function for /nl is plotted in Fig. 5. We 
found that the marginalized probability that /nl > was 
94%, the mean of /nl was 529 and the standard deviation 
194. 
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Figure 5. The marginalized posterior probability density func- 
tion function for /nl on cluster scales. 



We can combine the WMAP7 constraint / N l = 32 ±21 
(Komatsu et al. 2010) at wave-numbers of about 0.04 Mpc - 
to the one obtained here /nl = 529 ± 194 at wave-numbers 
of about 0.1 Mpc -1 to constrain the non-Gaussian spectral 
index, defined by 



r rCMB 

/nl = /nl 



fccMB 



(7) 



In order to use Eq. 7, we have to assume that /nl is the 
same sign on all scales. Using the CMB and cluster mean 
constraints on /nl implies ?ing ~ 3 which is quite high 
as it is generally first order in the slow roll parameters, but 
larger values are possible (Byrnes et al. 2009). Alternatively, 
such a rapid change in /nl could be explained by models 
that produce a step-function like change in /nl (Riotto and 
Sloth, 2010). 



4 CONCLUSIONS 

We found constraints of /nl = 529 ± 194 at wave-numbers 
of about 0.1 Mpc -1 using the observations of cluster 
XMM2235. We found p(/ NL > 0|XMM2235) « 0.94 and 
so XMM2235 is consistent with /nl = 0. The WMAP7 data 
was used to break the degeneracy with as- In order to do 
the analysis, we assumed XMM2235 was the most massive 
cluster of the survey. If this turns out not to be the case, our 
results will be conservative as a larger cluster mass would 
imply more preference for a higher value of /nl- 

As done in Jee et al. (2009), we based our analysis on 
the ACDM Jenkins et al. (2001) mass function fitted for the 
324 overdensity halo definition. We also tried several other 
mass functions. The mass of XMM2235 can be evaluated 
for other overdensities by assuming a NFW profile (Hu & 
Kravtsov 2003). We found similar results for the 180 over- 
density mass functions of Jenkins et al. (2001) and Sheth 
& Tormen (1999). It is a concern that in general the mass 
functions need to be extrapolated to some extent for the 
very high masses we are considering. This is partly due to 
the rareness of such massive halos in the /nl = case. It 
would be useful if this range could be mapped out better by 
a large number of low resolution N-body simulations. 

Our results are based on the assumption, that Matar- 
rese et al. (2000) approach for accounting for the enhance- 
ment of cluster numbers is valid up to /nl ~ 1000, which 
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should be tested against N-body simulations. Recent work 
(Enqvist et al. 2010) has shown that the high /nl cutoff in 
our marginalized distribution posterior distribution of /nl 
may be too drastic. This would imply that our results are 
conservative and that our mean /nl is higher. 

Also, selection effects could reduce the significance of 
our result in that rather than using the 11 deg 2 one should 
take the sky area of all surveys which have sufficient depth 
to see clusters at z ^ 1.4. For example if we take the 
the Hoyle et al (2010) compilation of high-z X-ray cluster 
surveys which covers an area of about 300 square degrees 
that would reduce our value of jj(/nl > 0|XMM2235) to 
about 15% although some of those surveys do not have good 
completeness for such a high redshift. Taking into account 
Sunyaev-Zeldovich surveys may reduce the significance fur- 
ther depending on whether they have sufficient depth to 
have seen clusters at z > 1.4. Additionally one could use 
a larger range of observed high-z clusters to constrain /nl. 
But a fuller analysis may require a detailed study using the 
selection functions of all current cluster surveys, so we leave 
this for future work. 

There are several other mechanisms, that could lead to 
unusually high mass clusters, such as features in the pri- 
mordial power spectrum (Chantavat et al. 2008), or modi- 
fications in the growth function due to dynamical dark en- 
ergy or non-Einstein gravity. Our results are consistent with 
/nl = but the allowed range of /nl is quite high. If our 
central values turn out to be correct, then future observa- 
tions should be able to detect the non-Gaussianity and its 
running (LoVerde et al. 2008, Jimenez and Verde 2009, Se- 
fusatti et al. 2009). 
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